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Abstract 

The self consistent version of the density functional theory is presented, which 
allows to calculate the ground state and dynamic properties of finite multi- 
electron systems. An exact functional equation for the effective interaction, 
from which one can construct the action functional, density functional, the 
response functions, and excitation spectra of the considered systems, is out- 
lined. In the context of the density functional theory we consider the single 
particle excitation spectra of electron systems and relate the single particle 
spectrum to the eigenvalues of the corresponding Kohn-Sham equations. We 
find that the single particle spectrum coincides neither with the eigenvalues 
of the Kohn-Sham equations nor with those of the Hartree-Fock equations. 
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I. INTRODUCTION 



The density functional theory (DFT), that originated from the pioneering work of Ho- 
henberg and Kohn [1], has been extremely effective in describing the ground state of finite 
many- electron systems. Such a success gave birth to many papers concerned with the gen- 
eralization of DFT, which would permit the description of the excitation spectra also. The 
generalization, on theoretical grounds, originated mainly from the Runge-Gross theorem, 
which helped to transform DFT into the time-dependent density functional theory TDDFT 
[3]. Both, DFT and TDDFT, are based on the one-to-one correspondence between particle 
densities of the considered systems and external potentials acting upon these particles. Un- 
fortunately, the one-to-one correspondence establishes only the existence of the functionals 
in principle, leaving aside a very important question on how one can construct them in re- 
ality. This is why the successes of DFT and TDDFT strongly depend upon the availability 
of good approximations for the functionals. This shortcoming was resolved to a large extent 
in [2,4,5] where exact equations connecting the action functional, effective interaction and 
linear response function were derived. But the linear response function, containing infor- 
mation of the particle-hole and collective excitations, does not directly present information 
about the single particle spectrum. 

In this Report, the self consistent version of the density functional theory is outlined, 
which allows to calculate the ground state and dynamic properties of finite multi-electron 
systems starting with the Coulomb interaction. An exact functional equation for the effective 
interaction, from which one can construct the action functional, density functional, the 
response functions, and excitation spectra of the considered systems, is presented. The 
effective interaction relating the linear response function of non-interacting particles to the 
exact linear response function is of finite radius and density dependent. We derive equations 
describing single particle excitations of multi-electron systems, using as a basis the exact 
functional equations, and show that single particle spectra do not coincide either with the 
eigenvalues of the Kohn-Sham equations or with those of the Hartree-Fock equations. 



II. EXACT EQUATION FOR THE FUNCTIONAL 

Let us briefly outline the equations for the exchange-correlation functional E xc [p] of the 
ground state energy and exchange-correlation functional Aj; C [p] of the action A[p] in the case 
when the system in question is not perturbed by an external field. In that case an equality 
holds [2,4] 

E xc [p] = A xc \p]\ p (r jU;=0 ), (1) 

since A xc is also defined in the static densities domain. The exchange-correlation functional 
E xc [p] is defined by the total energy functional E[p\ as 

E[p] = T k [p] + I [ P(ri)p(r f rfr^r 2 + E xc [p], (2) 

Z J 1 1*1 — T2 1 

where Tk[p] is the functional of the kinetic energy of the non-interacting Kohn-Sham par- 
ticles. The atomic system of units e = m = h= lis used in this paper. The exchange- 
correlation functional may be obtained from [2] 
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E xc [p] = -^J *2, iw, g') + 2np(r 1 )S(w)S(r 1 - r 2 )] 1^2 ^ ^ 

Equation (3) represents the expression for the exchange-correlation energy of a system [2] , 
expressed via the linear response function x( r i; r 2, iw, g'), with g' being the coupling con- 
stant. For Eq. (3) to describe A^fp] and -E^cfp] the only thing we need is the ability to 
calculate the functional derivatives of E xc [p] with respect to the density. According to Eq. 
(3), it means an ability to calculate the functional derivatives of the linear response function 
X with respect to the density p(r,u) which was developed in [2,5,6]. The linear response 
function is given by the integral equation 

X(r 1 ,r 2 ,u) = xo(ri,r 2 ,w) + J X o(r 1 ,r' 1 ,uj)R(r' 1 ,r' 2 ,uj)x(r' 2 ,r2,uj)dr' 1 dr' 2 , (4) 

with xo being the linear response function of non-interacting particles, moving in the single 
particle time- independent field [2,5]. It is evident that the linear response function x(d) 
tends to the linear response function of the system in question as g goes to 1. The exact 
functional equation for R(ri, r 2 , uj, g) is [2,5] 

R(ri, r 2 , uj, g) = - — r (5) 

1*1 — 1*2 



/ r x(r; ' 4 ™> 9>) w=^\ dri >dr2 ,d £ dg> - 



2 5p(r 1 ,u)5p(r 2 , -uj) 

Here R{r\, r 2 , uj, g) is the effective interaction depending on the coupling constant g of the 
Coulomb interaction. The coupling constant g in Eq. (5) is in the range (0 — 1). The single 
particle potential v xc , being time- independent, is determined by the relation [2,5], 

v xc(r) = J^ E *c[p\- (6) 

Here the functional derivative is calculated at p = p with p being the equilibrium density. 
By substituting (3) into (6), it can be shown that the single particle potential v xc has the 
proper asymptotic behavior [5,6], 

v xc (r -> 00) -> v x (r -> 00) -> -i. (7) 
The potential v xc determines the energies £j and the wave functions fa 

^~ + V H (r) + V ext (r) + v xc (r)j fa{v) = £i fa(r). (8) 

These constitute the linear response function Xo(*i, r 2, uj) entering Eq. (4) 

Xo = -n fc )0*(ri)0 i (r 2 )0fc(r 2 )0fc(ri) — — 1 

i,k 

and the real density of the system p, 



uj -uj ik + ir) uj + uj ik - ir] 



(9) 
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p(r)=5>|<Mr)| 2 . (10) 

i 

Here n» are the occupation numbers, V ex t contains all external single particle potentials of 
the system, viz. the Coulomb potentials of the nuclei. Eh is the Hartree energy 

2 J |ri-r 2 | 

with the Hartree potential V#(r) = SEh/Sp(t), and ou^ is the one-particle excitation energy 
^ik = £k — £i, and r] is an infinitesimally small positive number. 



III. THE EFFECTIVE INTERACTION 



The above equations (2-5) solve the problem of calculating E xc , the ground state energy 
and the particle-hole and collective excitation spectra of a system without resorting to 
approximations for E xc , based on additional and foreign inputs to the considered problem, 
such as found in calculations such as Monte Carlo simulations. We note, that using these 
approximations, one faces difficulties in constructing the effective interaction of finite radius 
and the linear response functions [1]. On the basis of the suggested approach, one can solve 
these problems. For instance, in the case of a homogeneous electron liquid it is possible to 
determine analytically an efficient approximate expression Rrpae for the effective interaction 
R, which essentially improves the well-known Random Phase Approximation by taking into 
account the exchange interaction of the electrons properly, thus forming the Random Phase 
Approximation with Exchange [4,5]. The corresponding expression for Rrpae is 



where 



RRPAE(q,g,p) = 



Re{q,9,P) = - 9 ~r 
Pf 



4ng 5E X 4ng 



+ 



\2p\ 



In 



Sp 



Ap 2 F 



+ R E (q,g,p), 



2pf 
3q 



In 



2p F - q 



2pi 



q 



l 

+ 3 



(12) 



(13) 



Here E x is the exchange energy given by Eq. (3) when x is replaced by xo- The electron 
density p is connected to the Fermi momentum by the ordinary relation p = p F /3n 2 . Having 
in hand the effective interaction in RRPAE^q, g, p), one can calculate the correlation energy e c 
per electron of an electron gas with the density r s . The dimensionless parameter r s = ro/as 
is usually introduced to characterize the density, with tq being the average distance between 
electrons, and ae is the Bohr radius. 



r s 


r c 
b M 


t RPA 


£ RPAE 


1 


-1.62 


-2.14 


-1.62 


3 


-1.01 


-1.44 


-1.02 


5 


-0.77 


-1.16 


-0.80 


10 


-0.51 


-0.84 


-0.56 


20 


-0.31 


-0.58 


-0.38 


50 


-0.16 


-0.35 


-0.22 
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In the Table, Monte Carlo results [7] e c M are compared with the results of the RPA calculation 
s c RPA , and e c RPAE when the effective interaction R was approximated by Rrpae [2,4]. The 
energies per electron are given in eV. Note that the effective interaction Rrpae{<1i p) permits 
the description of the electron gas correlation energy e c in an extremely broad range of the 
variation of the density. At r s = 10 the error is no more than 10% of the Monte Carlo 
calculations, while the result becomes almost exact at r s — 1 and is exact when r s — > [2,4]. 



IV. SINGLE-PARTICLE SPECTRUM 

Now let us calculate the single particle energies e*, that, generally speaking, do not 
coincide with the eigenvalues E{ of Eq. (8). Note that these eigenvalues Si do not have a 
physical meaning and cannot be regarded as the single-particle energies (see e.g. [1]). To 
calculate the single particle energies one can use the Landau equation [8] 

In order to illustrate how to calculate the single-particle energies q within the DFT, 
we choose the simplest case when the functional E xc is approximated by E x . As we shall 
see, the single-particle energies coincide neither with the eigenvalues calculated within the 
Hartree-Fock (HF) method nor with Si of Eq. (8). To proceed, we use a method developed in 
[5]. The linear response function xo and density p(r), given by Eqs. (9) and (10) respectively, 
depend upon the occupation numbers. Thus, one can consider the ground state energy E 
as a functional of the density and the occupation numbers 

E[p(r),ni] = T k \p(r),m] + E H [p{r),n t ] + E x \p(T),m] + J V ext (r)p(r)dr. (15) 

Here T k is the functional of the kinetic energy of noninteracting particles. As it follows from 
Eq. (3), the functional E x is given by [6] 

E x \p\ = ~ I Mr u r 2 ,iw) + 2np(r 1 )5(w)5(r 1 - r 2 )]^*^. (16) 

Upon using Eq. (16), the exact exchange potential v x (r) = 5E x /5p(r) of DFT can be 
calculated explicitly [6]. Substituting Eq. (15) into Eq. (14) and remembering that the 
single-particle wave functions fa and eigenvalues £j are given by Eq. (8) with v xc (r) = v x (r), 
we see that the single particle spectrum e$ can be represented by the expression 

o E x 

€i = Si- < fa\v x \fa > +- — . (17) 

orii 

The first and second terms on the right hand side in Eq. (17) are determined by the 
derivative of the functional T k with respect to the occupation numbers rij. To calculate the 
derivative we consider an auxiliary system of non- interacting particles in a field U(r). The 
ground state energy Eq of this system is given by the following equation 



EX 



T k + J U{v)p{v)dv. (18) 



5 



Varying Eq with respect to the occupation numbers, one gets the desired result 



Si = 



5El = 
Srij 5rij 



f < fa\u\fa >, 



(19) 



provided U — Vh + v x + V ext . The third term on the right hand side of Eq. (17) is related 
to the contribution coming from E x defined by Eq. (16). In the considered simplest case 
when we approximate the functional E xc by E x , the coupling constant g enters E x as a linear 
factor. If we omit the inter-electron interaction, g — > 0, that is, we put E x — > 0, we directly 
get from Eq. (17) — as it must be in the case of a noninteracting system of electrons. 
Note that it is not difficult to include the correlation energy in the simplest local density 
approximation 



E c [p,n,i\ = J p(r)e c (p(r))dr. 



Here the density p(r) is given by Eq. (10) and the correlation potential is defined as 

6E c [p] 



V c (r) = 



5p(r) 



(20) 



(21) 



Varying E[p(r), n«] with respect to the occupation numbers rii and after some straightforward 
calculations, we obtain the rather simple expression for the single particle spectrum 



e~ < fa\v x \fa > ~J2 n k J 



^r 1 ) ( /. i (r 2 )0*(r 2 )0 fc (r 1 ) 
ki - r 2 | 



dr 1 dr 2 . 



(22) 



Here £j are the eigenvalues of Eq. (8) with v xc = v x + V c . We employ Eq. (19) and choose the 
potential U as U — Vh + v x + V c + V ext to calculate the derivative 5T k /5rii. Approximating 
the correlation functional E c [p, n,] by Eq. (20), we simplify the calculations a lot, preserving 



at the same time the asymptotic condition, (y x + V c ] 



— 1/r. This condition is of 



crucial importance when calculating the wave functions and eigenvalues of vacant states 
within the framework of the DFT approach [5]. Note, that these functions and eigenvalues 
that enter Eq. (22) determine the single particle spectrum e^. This spectrum has to be 
compared with the experimental results. The single particle levels given by Eq. (22), 
resemble the eigenvalues ef F that are obtained within the HF approximation. If the wave 
functions <f>i would be solutions of the HF equations and the correlation potential V c (r) would 
be omitted, the energies €j would exactly coincide with the HF eigenvalues ef F . But this is 
not the case, since fa are solutions of Eq. (8), and the energies q do not coincide with either 
or with the eigenvalues Si of Eq. (8). We also anticipate that Eq. (22), when applied to 



HF 



calculations of many-electron systems such as atoms, clusters and molecules, will produce 
reasonable results for the energy gap separating the occupied and empty states. In the case 
of solids, we expect that the energy gap at various high-symmetry points in the Brillouin 
zone of semiconductors and dielectrics can also be reproduced. 



V. CONCLUSIONS 

We have presented the self consistent version of the density functional theory, which 
allows to calculate the ground state and dynamic properties of finite multi-electron systems. 
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An exact functional equation for the effective interaction, from which one can construct the 
action functional, density functional, the response functions and excitation spectra of the 
considered systems, has been outlined. We have shown that it is possible to calculate the 
single particle excitations within the framework of DFT. The developed equations permit 
the calculations of the single particle excitation spectra of any multielectron system such as 
atoms, molecules and clusters. We also anticipate also that these equations when applied to 
solids will produce quite reasonable results for the single particle spectra and energy gap at 
various high-symmetry points in the Brillouin zone of semiconductors and dielectrics. We 
have related the eigenvalues of the single particle Kohn-Sham equations to the real single 
particle spectrum. In the most straightforward case, when the exchange functional is treated 
rigorously while the correlation functional is taken in the local density approximation, the 
coupling equations are very simple. The single particle spectra do not coincide either with 
the eigenvalues of the Kohn-Sham equations or with those of the Hartree-Fock equations, 
even when the contribution coming from the correlation functional is omitted. 
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